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■, Abstract 

, We have studied the origin of the broad and skewed feature at 4.5-7.5 keV in the energy spectra of 

NGC 4151 using the ASCA and RXTE data. The feature consists of a narrow peak at 6.4 keV and a broad 
wing extended between 4.5-7.5 keV. An analysis of the long-term variations revealed that the feature 
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Introduction 



^ • became variable only on a time scale longer than 1.5 x 10^ s. Through a comparison with the continuum 

■ ' variabilities, wc found that the emission region of the excess flux at 4.5-7.5 kcV has an extent of 10^^ cm. 

The broad and skewed feature at 4.5-7.5 keV may be explained by the so-called "disk-line" model. If so, 
OO ' the size of the line-emitting region, 10^"^ cm, should be equal to several or ten-times the Schwarzschild 

radius of the central black hole. This results in a black hole mass of 10^^ Mq, which may be too large 
for NGC 4151. We propose an alternative explanation for the broad and skewed feature, i.e. a "reflection" 
model, which can also reproduce the overall energy spectra very well. In this model, cold matter with a 
• sufficiently large column density is irradiated by X-rays to produce a reflected continuum, which constitutes 

I the broad wing of the feature, and narrow fluorescent lines. The equivalent width of the iron fluorescent 

ff^ ' hne (~2 keV) and the upper limit of its width (cr < 92 eV) are also consistent with this model. From these 

' results and considerations, we conclude that the "disk-line" model has difficulty to explain the spectral 

^P , variations of NGC 4151, and the reflection model is more plausible. 

' Key words: galaxies: individual (NGC 4151) — galaxies: Seyfert — X-rays: galaxies — X-rays: 

individual (NGC 4151) — X-rays: sources 
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, Active Galactic Nuclei (AGN) are the most luminous objects in the universe. They emit large amounts of radiation 
over a wide range of wave-bands (from radio to 7-rays), and sometimes produce relativistic jets. It is a goal of modern 
astronomy to understand the origin of the extreme activities in AGN. 
' One of the key observational clues in the X-ray band to understand AGN is the iron line profile in the energy 
H I spectrum (Fabian et al. 2000). A broad and skewed line feature around 5-7 keV was clearly detected by ASCA for 
the first time in an energy spectrum of the Seyfert galaxy MCG-6-30-15 (Tanaka et al. 1995). Since then, a similar 
feature was detected by ASCA from several Seyfert galaxies (e.g. Nandra et al. 1997). The broad and skewed feature 
has been interpreted as an iron fiuorescent line originating from the innermost region of an accretion disk (Fabian et 
al. 1995, so-called "disk-line" model). In this model, a broad and skewed feature is explained by the combination of 
Doppler broadening due to the relativistic motion of the line emitting matter and gravitational redshift due to the 
central massive black hole. 

If the broad and skewed feature in the energy spectra is really a "disk-line" , we can determine the inclination 
angle of the accretion disk through an analysis of the line profile (Fabian et al. 1989). However, the thus-obtained 
inclination angle is sometimes inconsistent with that estimated with other methods (e.g. Nishiura et al. 1998). It may 
also be inconsistent with an expectation from the unified scheme, in which the accretion disk in Seyfert 2 galaxies 
should have an edge-on geometry. ASCA observations have revealed the presence of a broad and skewed line feature 
in several Seyfert 2 galaxies. The inclination angles have been deduced for these Seyfert galaxies from an analysis of 
the "disk- line", and always found to be around 30° (Turner et al. 1998). This inclination angle is unexpectedly small. 
This problem may be partially solved if we assume that the line is a composite and consists of a narrow line centered 
at 6.4 keV and a disk line with intermediate inclination angle (Weaver, Reynolds 1998). However, we still need an 
observational confirmation on such a composite model. 

NGC 4151 is a bright, nearby Seyfert galaxy (type 1.5; Osterbrock, Koski 1976; for a review see Ulrich 2000). The 
presence of the broad and skewed line feature in 4.5-7.5 keV has been known through ASCA observations (Yaqoob et 
al. 1995; Wang et al. 1999). Thus, NGC 4151 is one of the best targets to precisely study the line feature. A "disk-fine" 
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analysis of NGC 4151 has resulted in a face-on geometry of the accretion disk (Yaqoob et al. 1995), while observations 
of the [O III] A5007 image with the Hubble Space Telescope indicate an inclination angle of 60 ± 5° (Evans et al. 1993). 
This inconsistency might be partially resolved if we assume two disk lines and an additional narrow peak centered at 
6.4 kcV to reproduce the profile (Wang et al. 1999). However, it is not clear how the two disk lines with different 
inclination angles can be produced simultaneously. 

In this paper, we analyze the broad and skewed feature observed from NGC 4151 using a model-independent method 
as much as possible, and try to identify the origin of the feature. For this purpose, we analyze the long observations of 
NGC 4151 made with ASCA (in 2000) and RXTE (in 1999), while focusing on the spectral variations on various time 
scales. The results from another set of long ASCA observations of NGC 4151 have been presented by Wang, Wang, 
and Zhou (2001). 

2. Observations 

2.1. ASCA 

The primary ASCA data set used for the present analysis was acquired from 2000 May 13 through 25. Concerning the 
instrumentation on board ASCA, see the following references: Tanaka et al. (1994) for a general description of ASCA, 
Serlemitsos et al. (1995) for the X-ray telescope, Ohashi et al. (1996) and Makishima et al. (1996) for the Gas Imaging 
Spectrometer (CIS), and Burke et al. (1994) for the Solid-state Imaging Spectrometer (SIS). In the observations, SIS 
was operated in the 1-CCD Faint mode and CIS in the standard PH mode. A lower level discriminator was applied 
to SIS at 0.48 keV to avoid telemetry saturation due to flickering pixels. The data selection criteria for SIS and CIS 
are summarized in table 1. When the elevation angle, i.e. the angle between the field of view and the Earth limb, 
was less than 5°, the data were discarded because they were affected by scattered X-rays from the day Earth or by 
the absorption due to the atmosphere. Some of the optical light leakage is known in SIS, and the SIS data obtained 
when the elevation from the day Earth was less than 25° were discarded. The data with telemetry saturation are 
also discarded to avoid any deterioration of the detection efficiency. The net exposure of SIS was about 350 ks, which 
spans 1.1 X 10^ s, and an average count rate of SIS was 0.8 count s~^ after data screening and reduction. An X-ray 
source was detected with both SIS and CIS at a position consistent with the optical position of NGC 4151 within the 
positional accuracy of ASCA. A BL Lac object (1207-I-39W4; Arp 1997) was detected in the field of view of ASCA 
about 4.'7 north of NGC 4151. Contamination by this source was found to be negligible above 2 keV, which was used 
for the current data analysis. According to the Chandra observations, a faint source was also present at 2'.2 from 
NGC 4151 to the southwest (Yang et al. 2001). This source was not resolved with ASCA. Contamination by this 
source was completely negligible. 

The energy spectra and light curves of each detector were accumulated from a circular region centered on the 
source. The extraction radius of SIS was set to 3.'6, which is the maximum radius fit in both chips, on the whole. 
The background spectra of SIS arc accumulated from the source free region on the chip. The extraction radius of CIS 
was set to 6 arcmin, because various calibrations of CIS had been done with this standard extraction radius. The 
background data of CIS were accumulated from a circular region opposite to the target position against the bore-sight 
axis, excluding the source region. 

We usually obtained 4 sets of energy spectra from the ASCA data, two from SIS and the other two from CIS, for a 
single set of observations. To make the analysis simple, we summed up the two energy spectra from SIS and the two 
corresponding response matrices, respectively, which produced a single SIS spectrum with the corresponding response. 
Similarly, we summed up the two CIS spectra and the corresponding responses, respectively. The summed SIS and 
CIS spectra were fitted with a model spectrum simultaneously. The simultaneous fitting made a full utilization of 
the ASCA data in the sense that the higher detection efficiency below 5 keV and better energy resolution of SIS, and 
the higher detection efficiency above 5 keV of CIS are all reflected in the best-fit spectral parameters. Because there 
is 3% uncertainty of flux normalization between SIS and CIS, a systematic error of 3% was conservatively added to 
the SIS energy spectra in the simultaneous fitting. Although a decrement of SIS detection efficiency below 2 keV was 
reported by the instrument team (T. Yaqoob, private communications), we did not include the effect in the analysis 
because the spectral fittings were carried out only for the energy range above 2 keV. 

2.2. RXTE 

Data obtained by the Proportional Counter Array (PCA) on board the Rossi X-ray Timing Experiment (RXTE) were 
also used in the present analysis. Details of RXTE and PCA can be foimd in Jahoda et al. (1996). RXTE observations 
of NGC 4151 were carried out every '^5 days from the beginning of 1999 to study the long-term variations of the X-ray 
fiux. We used the data obtained from 1999 January 1 through July 24, which were publicly available from the archive 
in HEASARC^ at the time of data analysis. Reduction of the PCA data was carried out using the standard procedure 
(Edclson and Nandra 1999). In the analysis, we discarded data obtained when the elevation from the Earth limb was 

^ http://h6asarc.gsfc.nasa.gov/ 
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less than 10°, data obtained within 30 minute after the satelhte's passage of SAA, data obtained when the pointing 
offset was greater than 0.°02, and data when the count rate of electrons was greater than 0.1. The background spectra 
were calculated using a faint-source model according to a method described in the RXTE Cook Book.^ 

3. Analyses and Results 

3.1. ASCA Data 

Figure 1 shows the X-ray hght curve at 0.7-10 keV from NGC 4151 obtained with ASCA CIS on 2000 May 13-25. 
Here, the data are shown in 5.6 ks binning. This 5.6 ks bin-time corresponds to the orbital period of ASCA around 
the Earth. We chose this bin time because: (1) time variations shorter than the orbital period were not significant, 
(2) data gaps arc usually much shorter than the orbital period and we can get continuous light curve, (3) possible 
background variations, which are correlated with the orbital period, could be largely reduced. Note that data gaps 
were produced mainly by the earth occultation of the source, high background regions on the earth, and down-link to 
the ground stations. X-ray flux variations of a factor of 2-3 are clearly seen in the figure. 

We next studied the energy dependence of the flux variations. For this purpose, we calculated the ratios of the 
time-resolved energy spectra to the average energy spectrum. Considering the statistics of the data and the amplitude 
of the time variations, the time- resolved energy spectra were calculated every 1.8 x 10^ s and six spectra were obtained 
in total. We also calculated the average spectra for CIS and SIS, each of which correspond to the average of 6 
time-resolved spectra, as shown in figure 2. The ratios of the six time-resolved spectra to the average spectrum are 
presented in figure 3. Large time variations in both the flux and the shape can be seen at 2-5 keV, while the ratios 
are almost flat at 0.7-1.5 keV and 7-10 keV. The flat spectral ratios mean that the shape of the energy spectra does 
not vary at 0.7-1.5 keV and 7-10 keV. Furthermore, those at 0.7-1.5 keV are always close to unity. This indicates 
that the energy spectra change neither their shapes nor the fluxes in this energy range. It should be noted that a 
local structure is seen just around 6.4 keV and the ratios are closer to unity only in this narrow energy band. At 
least three spectral components may be required to reproduce these time variations of the spectral ratios: a stable 
soft component below ~2 keV, a variable hard component above ~2 keV, and a relatively stable component around 
6-7 keV. Because we are mostly interested in the feature around 6-7 keV, we do not discuss the soft stable component 
hereafter by limiting the energy range in the spectral analyses to above ^^2 keV. 

In order to study the nature of the relatively stable component around 6-7 keV, we need to separate the component 
from the underlying variable hard component. We first determined the shape of the variable hard component by 
performing a model fit to the average spectrum at 2-4 keV and 8-10 keV after masking the 4-8 keV range. The 
results show that a simple model of a power law with a single absorption cannot reproduce the spectrum, and that a 
"dual absorber" model is necessary. The need for a dual absorber for the ASCA spectrum was already pointed out by 
Weaver et al. (1994). The "dual absorber" model is expressed by the following formula: 

i^duai(-B) =A{CT.eM'<E)N^] + CT.)eM-<E)Nl]] E-^ , (1) 

where A is a normalization factor, F the photon index, A^h the column density, <j{E) the photo-electric absorp- 
tion cross-section, and Cp the source covering fraction (0 < Cp < 1). The superscript to A^h indicates two different 
column densities. In the following analysis, the photo-electric absorption cross-sections compiled by Morrison and 
McCammon (1983) are considered. The thus-determined continuum spectrum was subtracted from the average spec- 
trum; the residual feature at 4-8 keV is shown in figure 4. The presence of a narrow peak at 6.4 keV and a broad 
feature skewed toward the lower energy can be clearly seen. 

We then studied the time variations of the residual feature at around 6.4 keV. For this purpose, we used the 6 sets 
of SIS spectra, which have a better energy resolution than the CIS spectra. The continuum shape was determined 
separately for the 6 sets of spectra using the same method for the average spectrum. We set N"^ in equation (1) 
to to simplify the model fitting, considering the statistics of the data. We confirmed that setting N"^ to zero did 
not change the continuum shape at 4-8 keV more than 3% while taking the 4th data set (those with the largest 
fiux) as an example. The model was found to be acceptable for all 6 sets of the spectra with the largest value of 
X^/d.o.f. = 0.97/61. When the continuum model was interpolated to the 4-8 keV band, an excess feature could be 
clearly seen in all 6 sets of the energy spectra. We show them in figure 5 after subtracting the continuum model 
from the observed spectra. A prominent narrow peak at around 6.4 keV, and broad features at 4.5-6.0 keV (hereafter 
referred to as a red wing) and at 6.8-7.5 keV (a blue wing) are clearly seen in all 6 data sets. 

In order to study any profile change of the broad and skewed feature, we made the ratios of each of the 6 excess 
features to their average. The ratios were calculated by dividing the individual excess feature (after subtracting the 
interpolated continuum) by the excess feature in the average spectrum (figure 4). The results are shown in figure 6. 
By performing a test for the hypothesis that the ratios have no energy dependence, we found that the 6 sets of 
ratios are consistent with being constant in terms of the energy, and that the time variation of the profile of the broad 

^ http: //heasarc .gsf c .nasa. gov/docs/xte/recipes/cook_book.html 
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and skewed feature is not significant. At the same time, we calculated the flux ratio of the red wing (4.5-6.0 keV) to 
the narrow core (6.0-6.8 keV) and that of the blue wing (6.8-7.5 keV) to the narrow core for the 6 sets of SIS spectra. 
No significant time variation was found at the 90% confidence limit in the two sets of ratios (x^/d.o.f. = 1.22/5 at 
most). The standard deviation to the mean for the flux ratio of the red wing to the narrow core was calculated to be 
20% at most. Although the time variation is not significantly found, a variation of less than 20% in amplitude is not 
rejected. 

3.2. RXTE Data 

We analyzed the RXTE PCA data in order to study the spectral variations on a time scale longer than the 1.1 x 10® s 
covered by ASCA. We show in figure 7 the light curve calculated from the PCA data obtained from 1999 January 
1 through July 24. Significant flux variations can be clearly seen. The PCA spectrum averaged over the RXTE 
observations is shown in figure 8. As was done for the ASCA spectra, we fit a power law modified by "dual absorption" 
{N^ = 0) to the average spectrum in 2.8-10 kcV after masking the 4-8 kcV energy band. We also checked how it 
affected the continuum shape to set to zero, and found that there was no noticeable impact if the fit range was 
restricted to 2.8-4 and 8-10 keV. The ratio of the observed spectrum to the continuum model is shown in the lower 
panel of figure 8. We can clearly see the presence of a broad and skewed feature at 4-8 keV again. 

Next, in order to search for the time variations of the broad and skewed feature on time scales longer than 1.1 x 10® s, 
we divided the data into 11 subsets on a time bin of 1.5 x 10® s and calculated an energy spectrum for each subset 
of data. Wc then fit a power law modified by "dual absorption" to the 11 spectra in 2.8-10 kcV while excluding the 
4-8 keV band again. The fit was acceptable for all of the spectra. The residuals of the 11 spectra after subtracting the 
best- fit continuum model are shown in figure 9. A significant excess flux over the continuum model is always seen at 
4.5-7.5 keV. The excess has a peak at around 6.0-6.5 keV as well as a significant tail feature on the lower energy side. 
Then, we again investigated the time variations of the spectral shape of the excess by comparing the 11 excess spectra 
with their average. The spectral ratios of the 11 spectra to their average are shown in figure 10. We performed 
tests to check whether or not each of the 11 ratio-spectra is consistent with being flat. A flat model was found to be 
acceptable for all 11 ratios with x^/d.o.f. < 1.23/7, although the flux level of the excess shows a significant variation. 
This indicates that, although the excess flux at 4.5-7.5 keV is variable on a time scale longer than 1.5 x 10® s, the 
profile does not change significantly. 

4. Plausible Models 

In the previous section we discussed the time variability of the fiux and the profile of the broad and skewed feature 
consisting of a narrow peak at 6.4 keV and red/blue wings in a way fairly independent of any spectral models. In this 
section, we introduce two different model functions to reproduce the profile: one is the so-called "disk-line" model, 
and the other is a "reflection" model. 

4-.1. Disk-Line Model 

It is widely accepted that the broad and skewed feature at around 5-7 keV often seen from AGN is composed of 
gravitationally redshifted iron-fluorescent lines from matter orbiting in an accretion disk very close to the central black 
hole (Tanaka et al. 1995; Fabian et al. 1995). Hence, we first tried to fit the broad line-like feature observed from 
NGC 4151 with the disk-line model. We adopted a spectrum described by equation (1) as the continuum, where 
was set to to simplify the model. The outer and inner radii of the accretion disk were set to 1000 rsch and 10 rgch, 
respectively, where rsch is the Schwarzschild radius. The cmissivity of the line emission is assumed to have a power- law 
dependence on the radius with the index of —2 (Fabian et al. 1989). 

The model fitting was carried out in the energy range of 2.2-10 keV simultaneously to the SIS and GIS data for 
each of the 6 sets of the ASCA energy spectra. The photon index was fixed to 1.55, which was the average value when 
we fit a power law with a partial covering absorber to the spectra after masking 4-8 keV. The model was found to be 
acceptable for all of the spectra at 2.2-10.0 keV with x^/d.o.f. = (0.99— 1.10)/140. The best-fit parameters arc listed 
in table 2. The inclination angles and the line center energies obtained from the spectral fits to the 6 spectra arc not 
significantly different from one another. The reduced values for the hypothesis that the 6 values are constant were 
1.35 and 0.89 for 5 d.o.f. for the inclination angle and the line center energy, respectively. This means that these two 
parameters are consistent with being constant. The average inclination angle and the line center energy are 20° and 
6.44 keV, respectively. 

We applied the same model to the 11 energy spectra of PCA at 2.8-24 keV. However, the fit was not acceptable 
with x^/d.o.f. > 1.89/46. The deviation of the model from the observed spectra became noticeable above 16 kcV. 
Hence, we replaced the power law in equation (1) with a broken power law as described below: 
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where A is a normaUzation, Fi and F2 are photon indices, and i?bk is a break energy at which the spectral index 
changes. Note that the broken power-law is required due to the apparent change in the spectral slope above ^16 keV, 
not because we set A^^ to zero. 

The 11 energy spectra of PC A were all reproduced by the new model with xj/d.o.f. — (0.99 — 1.48)/46. The best- fit 
parameters are listed in table 3. We checked the time variation of the 11 values of each parameter by performing 
a test. The photon indices, Fi and F2, and the break energy, iSbki are found to have not changed significantly 
during the observational period. We cross-checked the spectral slope above the break energy, which was found to be 
F2 = 1.88 ±0.17, using the data from another detector HEXTE (High Energy X-ray Timing Experiment) on board 
RXTE. HEXTE consists of two clusters of 4 Nal/Csl phoswich scintillation detectors sensitive to X-rays from 15 
to 250 keV. The photon indices obtained from a simultaneous fit of a power-law model to the 11 sets of the PC A 
spectrum (18-24 kcV) and HEXTE spectrum (18-100 keV) are found to fall around 1.8 ±0.2. This is consistent with 
F2 1.9. It is noteworthy that, although the disk line model can always reproduce the excess feature around 5-7 keV, 
the continuum spectrum needs a break at around 17 keV to reproduce the energy spectrum up to 100 keV. 

^.2. Reflection Model 

In the disk-line model, the peak at 6.4 keV is considered to be a blue part of fluorescent lines coming from matter 
orbiting in an accretion disk very close to the central black hole. However, the peak energy is just the Kq, line of a 
neutral, or a lowly ionized iron, in the rest frame; furthermore, the profile around the peak can be reproduced by a 
single narrow line. Hence, it would be more natural to consider that the peak at 6.4 keV is a single fluorescent iron 
line coming from a region free from relativistic effects. If the line is emitted through the fluorescent process as a result 
of X-ray illumination on relatively cold matter far outside the central X-ray source. X-ray reflection at the surface 
of the matter should take place simultaneously with the fluorescent line emission. The X-rays reflected by the cold 
matter emitting the fluorescent lines could be a possible origin of the broad red and blue wings accompanied by the 
narrow line feature on both sides. Because X-ray reflection is due to Thomson scattering, the reflected X-rays should 
have the same spectrum as the illuminating X-rays in a typical X-ray band, but suffer from absorption when they pass 
through the cold mater. Thus, the reflected component may be approximated by applying cold-matter absorption to 
the continuum model. 

In order to test this possibility, we introduced the following model spectrum: 

F,,,{E) = P(^,F){C^exp[-a(^)iVi] + 1 - Ci^}exp[-a(^)iV2] + /(^ii„c) + C|P(^,F)exph(7(S)7V3]^ ^3) 

Here, the first term is the same as equation (1), representing continuum X-rays directly coming from the central X-ray 
source. The second term represents the fluorescent iron line at 6.4 keV and J(i?iinc) is a Gaussian function. The line 
width is assumed to be zero. The third term is for the reflection component. The efficiency of the X-ray reflection 
is represented by Cp, which corresponds to a covering fraction of the X-ray reflector optically thick for Thomson 
scattering. The function, P(i?, F), is a power law with a photon index of F. The photon index is optimized in the 
course of the fltting, because it may be affected by the introduction of a highly absorbed component. The column 
densities and covering fractions are also optimized. 

Model fitting was firstly carried out to the 11 PCA spectra in the energy range of 2.8-24 keV. In this energy range, 
because the effect of absorption by A^^ was negligible, we omitted this term. The fit was acceptable for all energy 
spectra with x^/d.o.f. = (0.87 — 0.99)/46. An example of the results of the model fitting is shown in figure 11 and the 
best-fit parameters are listed in table 4. The results show that 40-50% of the direct flux is reflected. We checked 
the constancy of the 11 values of each spectral parameter with a test. Wc found that the photon index and the 
absorption column of the reflected component, 7V^, arc both consistent with being constant (x^/d.o.f. = 0.98/10, and 
1.25/10, respectively). We conflrmed that the profile of the broad and skewed feature did not change significantly 
during the RXTE observations. 

The weighted mean of the 11 photon indices is F = 1.92±0.06. This is comparable to F = 1.8±0.2, which we obtained 
from a simultaneous fitting to the 11 sets of the PCA energy spectra at 18-24 keV and the HEXTE energy spectra at 
18-100 kcV. If we extrapolate the best-fit model at 2.8-24 keV to 100 keV, no significant discrepancy is recognized 
between the model and the observed energy spectra with HEXTE. This indicates that the intrinsic spectrum emitted 
from the central source should be a single power law with an index of ~1.9, and should extend up to 100 keV without 
a break, if the "reflection" model explains the broad and skewed feature at around 5-7 keV. Although the best-fit line 
energies are found to be consistent with being constant, the mean line center energy is i?iino = 6.1 ±0.1 keV. This is 
slightly lower than the center energy of the iron fluorescent line from neutral matter (6.4 keV). However, it is known 
that there is a systematic uncertainty of about 2-3% in the energy scale of the RXTE instrument. Thus, the line 
energy can still be interpreted as the iron fluorescent line. The time variation of the line flux is shown together with 
that of the continuum flux at 8-10 kcV in flgure 12. The time variation of the line flux is clearly seen in a time scale 
of 10^-10^ s. 

Next, we also applied the model function (equation 3) to the ASCA data. Because the energy range of the ASCA 
detectors is limited below 10 keV, some of the model parameters were not well constrained. Hence, we fixed the 
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following parameters in the fitting: F (= 1.9), iV^ {— 10^^ atom cm^^) and C| (= 47 %). Other parameters were 
optimized during the course of fitting. The model fitting was carried out to each of the 6 sets of the ASCA data, 
simultaneously to the SIS and GIS spectra. Fittings of this model were all acceptable with x^/d.o.f. = (0.94— 1.04)/140. 
An example of the results of the model fittings to the ASCA SIS spectra is shown in figure 13. The best-fit parameters 
are listed in table 5. The time histories of the line flux and the continuum flux at 8-10 keV for the ASCA data are 
plotted in figure 14. No significant change of the line flux can be seen on time scales of 10^-10^ s. As can be seen from 
the table, the values of N^, Cp and obtained from the spectral fits to the 6 sets of ASCA spectra significantly 
change on time scales of 10^-10^ s. The relative amplitude of the variations of and are ~20% and ^12%, 
respectively. 

Finally, we investigated the parameters of the narrow line component. We optimized the center energy and the 
intrinsic width of the narrow line component in the "reflection" model by fitting it to the average spectra of the ASCA 
observations in 2000. The line center energy and the intrinsic width were found to be 6.38 ±0.03 keV and less than 
92 eV (gaussian cr), respectively. 

4-3. Relation between the Line and the Continuum Fluxes 

An analysis of the time-resolved energy spectra showed that the iron line flux was consistent with being constant 
on time scales of 10^-10^ s, while significant variations were seen in its flux on time scales of 10^-10^ s. In order to 
study the origin of the line emission, the relation between the continuum flux and the line flux was studied. Here, 
the continuum flux was calculated in the energy band of 8-10 keV, which is higher than the K edge energy of neutral 
iron. In this energy range, the flux variation should reflect the variation of the intrinsic, power-law component. 

Because the line emission probably results from a reprocessing of the continuum X-rays by matter ambient to the 
continuum source, the light curve of the line flux might suffer from some amount of smearing and/or time delay to the 
continuum flux light curve. In fact, if we compare the light curves of the continuum flux and the line flux (figure 12), 
the line flux variation seems to follow the smeared variations of the continuum flux. In order to see the effect of 
smearing, we applied smearing to the continuum light curve using a following simple method. We redistributed the 
continuum flux of i-th time bin into N bins starting from i-th bin. Thus, the smeared light curve is expressed as 



where Fi is the flux of i-th bin in the original light curve, Gi the flux of i-th bin in the smeared light curve, and N 
the number of bins to smear. 

We then fitted the smeared light curve to the line flux light curve. We increased N from 1 until we obtain an 
acceptable fit. Note that = 1 corresponds to the case without smearing. The relation between N and x^/d.o.f. is 
plotted in figure 15. We could obtain an acceptable agreement with N = 4 between the smeared light curve of the 
continuum flux and the line flux light curve. The smeared light curve in the case of TV = 4 is compared with the light 
curve of the line flux in figure 16. The above-mentioned smearing algorithm of the continuum flux also introduces a 
delay of N/2 bins together with smearing of N/2 bins to both sides. Thus the typical delay and smearing time scale 
may be regarded as 3 x 10^ s. 

We then checked whether or not smearing of 6 x 10^ s can explain the absence of a correlation between the continuum 
flux and the line flux on shorter time scales for the ASCA data. Since the ASCA data is shorter than 6 x 10^ s, we 
could not smear the continuum flux light curve directly. Hence, we assumed that the observed light curve of the 
continuum flux in the ASCA observations repeats periodically. We then smeared the assumed light curve using the 
method described above. The fractional variation (i.e. standard deviation divided by the average) of the continuum 
flux was originally ~0.25, but was reduced to ~0.06 by smearing. The reduced value is consistent with the observed 
upper limit, 7%, of the fractional variation of the line flux in the ASCA observations. 

5. Discussion 

We analyzed the ASCA data obtained in 2000 May which covered time intervals of 1.1 x 10^ s. The data were 
divided into 6 sets with an integration time of 1.8 x 10^ s each to study the time variabilities of the spectrum. In order 
to analyze the excess component separately from the continuum, we masked a range of 4-8 kcV in spectral fitting and 
determined the continuum model. We found that a power law with a photon index of ~1.55 modified by two absorbers 
with different column densities and covering fractions can smoothly connect the two spectral parts at 2-4 keV and 
8-10 keV. Above this continuum, an excess is clearly detected in 4.5-7.5 keV, which has a broad and skewed feature. 
The feature has a prominent narrow peak at 6.4 keV, but accompanies a broad red wing at 4.5-6.0 keV and a blue 
wing at 6.8-7.5 keV. The flux and shape of the excess feature obtained from each of the 6 spectra were compared 
with one another. It is found that the excess flux is consistent with being stable, concerning both the flux and the 
spectral shape, over the observations of 1.1 x 10^ s. 




(4) 



j=i-N+l 
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We also analyzed data obtained by RXTE from January through July in 1999 in order to search for any time 
variations of the iron line on a time scale of > 10^ s. We obtained 11 sets of time-resolved energy spectra, each of 
which covers a time interval of 1.5 x 10^ s. If we selected only the energy bands free from the iron structures, i.e. the 
2.8-4.0 kcV and 8-10 kcV bands, each energy spectrum could be reproduced again by a power law with a photon index 
of ^1.5 modified by the partial covering absorption. This result is consistent with that of ASCA. In 4.5-7.5 keV, 
an excess flux over the continuum is clearly seen as a broad feature, in which significant flux variations are detected. 
However, no clear change was noticed in its profile. The absence of a profile change in spite of the significant flux 
variations in the broad and skewed feature strongly indicates that the feature is produced through a single mechanism. 

We studied the time scale of the flux variations of the narrow line. It was found that the flux variations are not 
significant on time scales of 10^-10^ s, whereas they become significant on time scales of 10^-10^ s. We also studied 
whether or not an introduction of a smearing effect can improve the correlation between the line flux and the continuum 
flux. Wc redistributed the continuum flux light curve with a simple method approximating a smear and delay. We 
found that the agreement between the line flux light curve and the smeared light curve of the continuum flux becomes 
acceptable when we introduce a smear and a delay on a time scale of 3 x 10^ s. These results concerning the line flux 
variations and the effect of a smear strongly suggest that the line emitting region should have a size extent of 10^^ cm. 

We further studied the time scale of the variation in the absorption column density. A column density of the order 
of ~10^^ atom cm~^ is necessary to reproduce the continuum spectrum, and is found to vary significantly on a time 
scale of 10^-10^ s with a relative amplitude of ~20%. Taking account of the relative amplitude of the variation, the 
size of the absorber is indicated to be no larger than 10^^ cm, which is just the size of the line emitting region. The 
absorber should be around, or inside, the line emitting region. 

We tried a continuum plus a disk-line model to reproduce the energy spectra, and found that this model can 
reproduce the energy spectra at 2-10 keV for both ASCA and RXTE data. As a result of the previous discussion, 
the size of the line emitting region should be as large as 10^^ cm. In the disk-line model, the line is assumed to be 
emitted from a region with a size of several to ten times the Schwarzschild radius. If the size of 10^^ cm deduced from 
observational results corresponds to several to ten-times the Schwarzschild radius, wc need to assume a central mass 
close to 10^^ Mq. This central mass is not consistent with 10^ Mq estimated for NGC 4151 from various methods, 
e.g. X-ray variability time scale (Hayashida et al. 1998) and the analysis of the high ionization lines in the broad line 
region (Clavel et al. 1987; Ulrich, Horne 1996). This mass is rather close to the typical mass of a galaxy, and would 
be too large as the mass of the central black hole. 

The presence of the absorber around or inside the line emitting region would also be difficult to explain in the 
context of the "disk-line" model. If the line emitting region has a size of several to ten-times the Schwarzschild radius, 
the absorber should necessarily be located at a region very close to the central black hole. This may not be consistent 
with the Unified scheme in which the heavy absorption as seen in NGC 4151 is considered to be due to a dusty torus 
around the central active region. We obtain an inclination angle of the disk to be ~20°. This inclination angle is 
smaller than the generally accepted value 63° , which is supported by the Chandra observation (Ogle et al. 2000) and 
the optical observations. This is also a disadvantage of the "disk-line" model. 

As discussed above, there exist some serious difficulties in the disk-line model, which are not consistent with a 
reasonable, common picture of AGN, although we cannot completely exclude its possibility. We consider here an 
alternative explanation for the origin of the broad and skewed profile, in which the narrow peak at 6.4 keV should be a 
narrow line and the red and blue wings should be a part of the continuum. However, according to the constant profile 
of the broad feature around 5-7 keV, the narrow line and the continuum should have a strong physical connection. 
The line energy at 6.4 keV strongly suggests that the line is the fluorescent K-linc from neutral or low ionization 
iron. If the line is really emitted through a fluorescent process, it implies that some amount of matter exists in the 
vicinity of the X-ray source. This matter should be irradiated by X-rays from the X-ray source and the fluorescent 
lines should be emitted there. If this is the case, continuum X-rays should also be emitted from the matter through 
electron scattering of irradiating X-rays. This reflected component should be observed together with the fluorescent 
line, and could be the red and blue wings. 

When an X-ray from the source is absorbed by an iron atom in the X-ray reflector through photo-ionization of a 
K-electron, a fluorescent K-line is emitted with a certain probability called a fluorescent yield. On the other hand, 
when an X-ray from the source hits an electron in the reflector, it could be re-cmitted towards us through Thomson 
(Compton) scattering. The cross section of Thomson scattering is about 10^^^ cm^. Hence, if the column density of 
the X-ray reflector is sufflciently large, an X-ray penetrated into the reflector would, on average, experiences photo- 
electric absorption by matter with a column density of about 10^^ cm~^ before being scattered by an electron. As a 
result, X-rays reflected by sufficiently thick matter should have a spectrum with a continuum shape being the same as 
the intrinsic X-rays, but suffering from photo-electric absorption by matter with a column density of about 10^^ cm~^. 
In fact, the broad and skewed feature can be reproduced by a model consisting of a narrow line at 6.4 keV, and 
a power law with the same slope as the remaining continuum and with a photo-electric absorption by a column of 
10^^ cm~^. The equivalent width of the 6.4 kcV line with respect to the heavily absorbed component is about 2 keV. 
This is roughly consistent with a value expected from a case when the heavily absorbed component is just the reflected 



8 



K. Takahashi, H. Inoue, and T. Dotani 



component (see e.g. Inoue 1985). 

If we adopt the above model for the broad and skewed feature around 4.5-8 keV, the total spectrum becomes to have 
three differently absorbed components with the same continuum shape. The absorption column densities are about 
5 X 10^^ cm~^, (1-2) X 10^'^ cm~^, and about 10^^ cm~^. Among the three components, two components absorbed by 
smaller columns are found to be variable on a time scale of 10^-10^ s, whereas the component with the absorption 
by ~ 10^'^ cm~^ was steady on that time scale. This strongly suggests that the heavily absorbed component has a 
different origin from the other two components. The time variability of the heavily absorbed component is rather 
similar to that of the fluorescent iron line. This indicates a strong physical coupling between the fluorescent iron 
line and the heavily absorbed component, and strongly supports the idea that the heavily absorbed component is a 
reflected emission by cold matter which is also responsible for the fluorescent iron line. 

The above model fits both the ASCA spectra at 2-10 keV and the RXTE spectra at 2-24 keV weU. The model 
can also fit the data in a higher energy range. Actually, if we extrapolate the best-fit model in the energy range of 
2.8-24 keV to 100 keV, no significant discrepancy is recognized between the model and the observed energy spectra. 
Therefore, this model can explain wider range of energy spectra without the spectral break than the model based 
on the disk-line. If the reflection is really at work, the covering fraction of the absorber may be related to the solid 
angle of the reflector, 17, subtending to the X-ray source. The covering fraction that we obtained for an absorber of 
10^^ cm~^ is '^0.47, which corresponds to ~ 0.47 x Air. The reflection structure was detected in the energy spectra 
of NGC 4151 with the Ginga/OSSE data (Zdziarski et al. 1996) and the BeppoSAX data (Piro et al. 1998). The 
BeppoSAX data indicated that the relative contribution of the reflection changed with time. The covering fraction 
of the reflector so far reported ranges over fl ^ (0.0 — 0.4) x in (Zdziarski et al. 1996; Piro et al. 1998). The covering 
fraction which we obtained is slightly larger than this range, but is not very different, if we consider the different 
modeling of the reflection component. 

Based on the reflection model, we may be able to constrain the mass of the central black hole using the iron line 
parameters. From the model fitting, we found that the iron line is consistent with having no intrinsic width, and its 
upper limit is 92 eV. We also found that the iron line emitting region has a size of 10^^ cm. It may be natural to 
assume that the line emission region is located around the central black hole and is rotating at the Kepler velocity 
around the black hole. Then, the line width may be determined by the Doppler effect; the upper limit of the line 
width can be related to the upper limit of the Kepler velocity as w/c ~ 0.014. This may be converted to the mass of 
the black hole using the relation ~ GM/r, which yields the upper limit of the black hole mass, ~10^ Mq. Here, 
we assume r 10^'^ cm. This is consistent with the estimated mass of 10^ Mq. 

If the mass of the central black hole is 10*" Mq, a few times 10^^ cm corresponds to 10^ fsch- At this distance, the 
X-ray reflector should exist. Since the reflected component shows evidence of absorption by a column of lO^"' cm~^, 
the reflector should be Thomson thick. The solid angle of the reflector, as viewed from the central X-ray source, 
should be about 0.47 x in. This reflector would be a so-called dust torus, which is generally believed to exist around 
the central engine in Seyfert galaxies. 

It has been shown that the continuum spectrum of NGC 4151 needs a partial covering by X-ray absorbers with a 
column density of (1-2) x 10^'^ cm^^. The absorption column and the covering factor have been found to vary on a 
time scale of 10^-10^ s. As already discussed above, this leaky absorber should exist around or inside the line emitting 
region. If a dusty torus emits the line, a geometrical relation between a dusty torus and a broad line region seems to 
be consistent with the observed constraint on the geometrical relation between the line emitting region and the leaky 
absorber. The absorbers could be relatively cold clouds in the broad line region. If so, the typical time scale of a 
change of the leaky absorption should be roughly estimated by R/V, where R is the size of the X-ray emitting region 
and V the velocity of the X-ray absorbing clouds. Since the typical velocity of the broad line clouds is several times 
10^ cm s~^ and the typical time scale of the absorption change is a few times 10^ s, the size of the emission region 
may be about 10^^ cm. This size corresponds to a few ten-times rsch of a black hole with a mass of 10^ Mq, and is 
consistent with a natural expectation that a region close to the central black hole should be emitting the continuum 
X-rays. 

From these results and considerations, wc conclude that the disk line model has difficulty to explain the spectral 
variations of NGC 4151, and the reflection model is much more plausible. 
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Table 1. Selection criteria of ASCA data. 



Telemetry bit rate High and Medium 

Pointing offset from the target < 0.05° 

Elevation angle from the Earth limb > 5° 

Cut-off rigidity > 6 GV 

Elevation angle from the day Earth (SIS only) > 25° 



Table 2. Best-fit parameters of the "disk-line" model for the 6 energy spectra from ASCA. 



Sequence number 


1 


2 


3 


4 


5 


6 


[W^ atom cm~^] 


11±2 


9±3 


7±2 


5±3 


7±3 


7±3 


cM%] 


87±4 


90 ±4 


90 ±4 


93±4 


87 ±5 


84 ±5 


Inclination angle [deg] 


21±5 


18 ±4 


21 ±5 


22 ±6 


21 ±6 


19 ±6 


line flux [10~^ photon 


cm-2 s-i] 3.0 ±0.3 


3.0 ±0.3 


3.0 ±0.3 


3.5 ±0.3 


2.8 ±0.3 


3.0 ±0.3 


Xl/d.ol 


1.10/140 


1.02/140 


0.99/140 


1.03/140 


1.02/140 


1.03/140 


Note: See equation (1) 


1 for the continuum model. Errors are at the 90% confidei 


Qce limit. 





Table 3. Best-fit parameters of the "disk-line" model for the 11 energy spectra from RXTE PCA. 



Sequence number 


1 


2 


3 


4 


5 


6 


Photon index (Fi) 


1.50±0.11 


1.45 ±0.13 


1.58±0.11 


1.60±0.11 


^•^^-0.12 


1 4^+0.10 

^•^0-0.12 


Break energy [keV] 


17±2 


18 ±2 


17±2 


18 ±2 


16 ±3 


19 ±3 


Photon index (r2) 


1.92 ±0.22 


^•"^-0.18 


1.89 ±0.21 


2.01 ±0.23 


1.0»_g 22 


1 87+0.24 
^•^''^-0.26 


A^^ [10^2 atom cm^^j 






4.9 ±2.2 


6 1+2° 


3.4±2.1 


3.8±2.2 





82 ±15 


> 70 


83 ±17 


75 ±18 


>68 


> 79 


Line flux [10 photon cm ^ s ^] 


8.9 ±0.6 


9.1±0.6 


9.8 ±0.6 


9.4±0.6 


9.0 ±0.6 


8.0±0.6 


X?/d.o.f 


1.07/46 


1.40/46 


1.12/46 


1.48/46 


1.14/46 


1.27/46 


Sequence number 


7 


8 


9 


10 


11 




Photon index (Pi) 


1.51±0.12 


1 t;i+0.ib 
J^-"JJ^-0.13 


1.50±0.14 


1.45 ±0.13 


1.45±0.11 




Break energy [keV] 

Photon index (r2) 

[10^2 atom cm^^] 

c-i [%] 

Line flux [10^^ photon cm^^ s^^] 


17±2 
1 q2+°-2i 

^■^^-0.24 

5.3 ±2.3 
87 ±12 
7.5 ±0.6 


16 ±2 

^•^''-0.28 

7 4+2-5 

'•^-2.7 

6.5±0.6 


17±3 

1.0»_Q 24 

6.0 ±0.6 


19 ±2 

7 7+2.0 
'•'-2.4 

6.7±0.6 


16 ±2 

1 00+0.26 
i.»^_0.24 

4.5±2.3 

> 78 
6.1±0.5 




X^/d.o.f 


1.12/46 


0.99/46 


1.39/46 


1.28/46 


1.20/46 





Note: See equation (1) for the continuum model. Errors are at the 90% confidence limit. 
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Table 4. Best-fit parameters of the "reflection" model for the 11 energy spectra from RXTE PC A. 



Sequence number 1 2 3 4 5 6 

Photon index 1.90±0.12 1.91+^^3 1.97±0.10 ^.OOtois 1.87±0.15 1.89±0.10 

iV^ [10^2 atom cm-2] 12 ±2 11 ±2 10 ± 3 14 ±2 10 ± 3 12 ±2 

[%] 84±6 82±4 77±6 77±4 78±8 75±7 

A^i^ [10^2 atom cm-2] IO6+2? 109 ± 27 99 +22 ^^6^11 90 ±30 117l2? 

[%] 53±13 56±12 52±11 61 ± 13 53±18 37±10 

Line flux [10--* photon cm-2 s-i] 4.9 ±0.6 5.1 ±0.6 5.7 ±0.7 5.4 ±0.6 5.3 ±0.7 4.3 ±0.7 

xg/d.o.f 0.99/46 0.94/46 0.91/46 0.88/46 0.97/46 0.87/46 

Sequence number 7 8 9 10 11 

Photon index 1.89+j;-f2 1.82+'^-}!^ 1.85 ±0.14 1.89±0.13 1.83 ±0.12 

iV^ [10^2 atom cm-2] 12 ±1 14 ±2 13 ± 2 18 ±2 12 ± 1 

[%] 76±5 81±6 83±8 85±4 82±6 

iV| [1022 atom cm-2] 120 ± 34 110 ±35 90 ± 33 100 ± 34 lOOtg? 

C| [%] 35 ±10 33 ±11 43^19 66 ±13 45 ±11 

Line flux [10^-^ photon cm-2 s-i] 4.4 ±0.6 3.8 ±0.6 3.8 ±0.6 3.9 ±0.6 3.2 ±0.5 

xl/d.o.f 0.98/46 0.91/46 0.99/46 0.91/46 0.90/46 



Note: Errors are at the 90% confidence hmit. 



Table 5. Best-fit parameters of the "refiection" model for the 6 energy spectra from ASCA. 



Sequence number 1 2 3 4 5 6 

[1022 atom cm-2] r7±2 15lb3 16±3 12±3 18±3 19 ±3 

C^[%] 67±5 59±7 50±7 46±5 47±5 48±5 

iVg [1022 atom cm-2] 4.6 ±0.4 5.3 ±0.5 5.6 ±0.5 4.9 ±0.6 4.9 ±0.5 4.6 ±0.5 

fine flux [lO"'' photon cm-2 s-i] 1.9 ±0.2 1.9 ±0.3 1.8 ±0.3 1.9 ±0.3 1.8 ±0.3 2.0 ±0.3 

xg/d.o.f 1.04/140 0.97/140 0.94/140 0.99/140 0.99/140 1.02/140 

Note: The model of equation (3) is used. The photon index, the column density of iV^, and the 
covering fraction of C| are fixed to 1.9, 102'^ atom cm^2^ a^d 47%, respectively. Errors are at the 
90% confidence limit. 
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Fig. 1. X-ray light curve of NGC 4151 in 5.6 ks time bin calculated with the ASCA GIS data obtained during 2000 May 13-25. 
The X-ray flux varied by a factor of 2-3 during the observations. 
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Fig. 2. Average energy spectra (folded with the detector response) of NGC 4151 obtained with ASCA in 2000 May. The spectra 
of GIS (crosses with open circles) and SIS (bare crosses) are shown separately. A prominent line like structure can be seen around 
6.4 keV in both energy spectra. 
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Fig. 3. Ratios of the time-resolved energy spectra (obtained every 1.8 X 10^ s) to the average energy spectrum. The left panel shows 
those obtained from GIS and the right panel from SIS. The ratios are consistent between GIS and SIS. Time sequence of the spectra 
is indicated in the left panel. 
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Fig. 4. Excess feature in 4-8 keV after subtracting the continuum model. The continuum model was determined by fitting the 
average spectrum in two energy ranges, 2.2—4.0 keV and 8-10 keV. The data with open circles and those with filled squares show 
the results for SIS and GIS, respectively. 



14 



K. Takahashi, H. Inouc, and T. Dotani 



[Vol. , 



I 

> 

o 
U 



0.02 
0.01 


0.02 
0.01 


0.02 
0.01 




t ' 

-++d^++ !^VVt*.^!^-.. 


t ' ' : 
- ^++++"'* \ ' 

-H +^+^ 
-+^rt_^+ +t>W^rvt^v - 


' 1 ■ 1 ' 1 ' 1 — 

i T 

1+++++^+ ^.yiV-^«'*^-^- - 


— ' — 1 — ' — 1 — ' — 1 — ' — 1 — 

5 

- f+-H-[ -|- . t+<Wfn+ti.H,+. 




1 +Hrt^4^T %t^+tV; W . . 

■ 1 I 1 ■ 1 ■ 1 


6 It 



8 10 4 6 

Energy (keV) 



8 10 



Fig. 5. Profiles of tlie excess features in 4-8 keV of the six time-resolved SIS spectra after subtracting the continuum model. The 
time sequence is indicated in the panels. 
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Fig. 6. Ratios of the excess features (see figure 5) to their average deduced from the ASCA SIS data. All of the ratios are consistent 
to be constant. This means that the excess feature in 4-8 keV show little time variations. 
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Fig. 7. Light curve of NGC 4151 in 8-24 keV calculated from the RXTE PGA data obtained from 1999 January 1 through July 
24. Significant flux variations are clearly seen. 




Fig. 8. PGA spectrum of NGC 4151 averaged over the RXTE observations from 1999 January 1 through July 24, (crosses in the 
upper panel) together with the best fit continuum model to the spectrum after masking 4—8 keV range (histogram in the upper 
panel). The lower panel shows the ratio of the observed spectrum to the continuum model. 
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Fig. 9. Residuals of the 11 PCA spectra after subtracting the continuum models. The residuals were calculated from spectra 
accumulated for 1.5 X 10® s each by subtracting the continuum model determined by excluding the 4-8 keV band. A significant 
excess flux over the continuum model is always seen in the 4.5-7.5 keV band. 




Fig. 10. Ratios of the 11 residual spectra (see figure 9) to their average. The ratios are all consistent with being constant, although 
the average flux level (indicated by the broken line) changes among the 11 ratios. 



Origin of the "Disk-Line" Speetrum Feature in NGC 4151 



17 




Energy (keV) 



Fig. 11. Example of the fits of the "reflection" model to the PCA speetra. The erosses represent the observed speetrum, while the 
histograms represent the two model-components, the direct power-law component and the reflected power law -|- the narrow line 
component. 
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Fig. 12. Light curves of the continuum flux (upper panel) and the iron line flux (lower panel) on a 1.5 X 10^ s time bin. These light 
curves were obtained form RXTE observations. The continuum flux was calculated in 8.0—10 keV. 
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Fig. 13. Spectral fit of the "reflection" model to the ASCA SIS spectrum of NGC 4151. The crosses represent the observed 
spectrum. The histograms are two model-components: the direct power-law component and the reflected power law -|- the narrow 
line component. The lower panel shows the ratio of the observed spectrum to the model-spectrum. 
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Fig. 14. Light curves of the line flux (lower panel) and the continuum flux (upper panel) with 1 cr error are plotted in 1.8 X 10^ s 
binning. These curves are obtained from the ASCA data. The continuum flux was calculated at 8.0-10 keV. The line flux docs not 
show significant time variations. 
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Fig. 15. Reduced value when we compare the smeared light curve of the continuum flux with that of the line flux, as a function 
of a parameter, N , representing the degree of smearing. The reduced value becomes acceptable when A'^ increases up to 4. 
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Fig. 16. Light curves of the line flux and the smeared continuum flux. For the calculation method of smearing, sec text. A good 
agreement can be seen between them. 



